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Abstract 

We consider the computation of the mean of sequences in the 
quantum model of computation. We determine the query complex- 
ity in the case of sequences which satisfy a p-summability condition 
for 1 < p < 2. This settles a problem left open in Heinrich (2001). 

1 Introduction 

Computation of the mean of sequences and, equivalently, summation of se- 
quences, is an important numerical task, in particular for huge number of 
summands occurring in many numerical applications such as, e.g., high di- 
mensional integration. The larger the number of summands (the larger the 
dimension), the less these problems are tractable on a classical computer. 
It is therefore an interesting and challenging task to understand to which 
extent a quantum computer could bring speed-ups. First results for the 
summation of bounded sequences are due to Grover (1998), Nayak and Wu 
(1999), Brassard, H0yer, Mosca, and Tapp (2000). The case of sequences 
satisfying a p-summability condition, which arises in various problems like 
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integration of functions from Lp and Sobolev classes, was studied in Hein- 
rich (2001). Up to logarithmic factors for p = 2, in the case 2 < p < oo the 
query complexity of tlic summation problem was determined. For the case 
1 < p < 2, matching upper and lower bounds were obtained only under an 
additional restriction. The bounds for the remaining case did not match. 
In this paper we settle this problem and determine the query complexity in 
the full range of parameters. 

Applications of our results to the quantum complexity of integration of 
functions from Sobolev classes are given in Heinrich (2001a). The use of 
quantum summation for integration was first pointed out by Abrams and 
Williams (1999). The quantum complexity of integration was studied in 
Novak (2001), later in Heinrich (2001) and Heinrich and Novak (2001). Path 
integration is discussed in Traub and Wozniakowski (2001). Furthermore, 
we refer to the surveys Ekert, Hayden, and Inamori (2000), Shor (2000), 
and to the monographs Pittenger (1999), Gruska (1999) and Nielsen and 
Chuang (2000) for general reading on quantum computation. 

Our analysis is based on the framework introduced in Heinrich (2001) 
of quantum algorithms for the approximate solution of problems of analy- 
sis. This approach is an extension of the framework of information-based 
complexity theory (see Traub, Wasilkowski, and Wozniakowski, 1988, No- 
vak, 1988, and, more formally, Novak, 1995) to quantum computation. It 
also extends the binary black box model of quantum computation (see, e.g., 
Deals, Buhrman, Cleve, and Mosca, 1998) to situations where mappings 
from spaces of functions to the scalar field (such as the mean or the inte- 
gral) have to be computed. Let us recall the main notions here. For more 
details and background discussion we refer to Heinrich (2001). 

Let D, K be nonempty sets, let J^{D,K) denote the set of all functions 
from D to K, and let F C J-'{D, K) be a nonempty subset. Let K, the 
scalar field, be either R or C, the field of real or complex numbers, let G 
be a normed space over K, and let S" : F — ^ G be a mapping. We seek to 
approximate S{f ) for / € F by means of quantum computations. Let Hi 
be the two-dimensional complex Hilbert space C^, with its unit vector basis 
{eo,ei}, let 

be the tensor product of m copies of Hi, endowed with the tensor Hilbert 
space structure. The following notation is convenient: 

Z[0,Ar) := {0,...,iV-l} 

for iV G N (as usual, N = {1,2, . . . } and Nq = Nu{0}). Let Cm = : i e 
Z[0, 2"^)} be the canonical basis of iJ^, where |i) stands for e^Q (8) • • • (8) , 
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^ — YlT=o 3k'^^~^~^ the binary expansion of i. Denote the set of unitary 
operators on Hjn by U{Hm)- 

A quantum query on F is given by a tuple 

Q = {m,m',m",Z,T,l3), (1) 

where m,m',m" G N, m' + m" < m,Z C Z[0, 2™') is a nonempty subset, 
and 

t: Z D 

are arbitrary mappings. Denote m{Q) := m, the number of qubits of Q. 

Given such a query Q, we define for each f G F the unitary operator Qf 
by setting for \i) \x)\y) £ Cm = Cm' ® Cm" ® Cn 



-"m—m' ~m" 



<V/ \^) \x) W - <y I otherwise, ^ ' 

where © means addition modulo 2™". 

A quantum algorithm on F with no measurement is a tuple 

A = (Q,(i7,)^=o), 

where Q is a quantum query on F, n G No and Uj G U{Hm) (j = 0, . . . , n), 
with m = m{Q). Given / G F, we let Af G U{Hm) be defined as 

Af = UnQfUn-l...UiQfUo. (3) 

We denote by nq{A) := n the number of queries and by m{A) = m = m{Q) 
the number of qubits of A. Let (^/(a^, y))x,yeZ[o,2"») be the matrix of the 
transformation ^/ in the canonical basis Cm- 

A quantum algorithm on F with output in G (or shortly, from F to G) 
with A; measurements is a tuple 

A = {{A,)1-',,{h)1-l^l 

where A: G N, and A(_{1 = {), . . . ,k — 1) are quantum algorithms on F with 
no measurements, 

6oGZ[0,2"*°), 
for 1 < £ < A; — 1 , 6^ is a function 

6,:nZ[0,2"^^)^Z[0,2-0, 

j=0 
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where we denoted := m{A(), and </? is a function with values in G 

k-l 

e=o 

The output of A at input f € F will be a probability measure A{f) on G, 
defined as follows: First put 

PA,f{xo,---,Xk-i) = \Aoj{xo,bo)f\Aij{xi,bi{xo))\'^ ... 

■ ■ ■ \Ak-i,f{xk-i, bk-i{xo, Xk-2))f. (4) 

Then define A{f) by setting for any subset C ^ G 

Mf){C) = Yl PAj{xQ, Xk-i). (5) 

ip{xo,-,Xk-i)&C 

By nq{A) := ^^^q nq{A(:) we denote the number of queries used by A. 

Informally, such an algorithm A starts with a fixed basis state 60 and, 
at input /, applies in an alternating way unitary transformations Uoj (not 
depending on /) and the operator Qf of a, certain query. After a fixed 
number of steps the resulting state is measured, which gives a (random) 
basis state, say ^o- This state is memorized and then transformed (e.g., by 
a classical computation, which is symbolized by 61) into a new basis state 
61 (^o)- This is the starting state to which the next sequence of quantum 
operations is applied (with possibly another query and number of qubits). 
The resulting state is again measured, which gives the (random) basis state 
^1. This state is memorized, 62(^0) ■^i) is computed (classically), and so 
on. After k such cycles, we obtain • • • iCfc-i- Then finally an element 
of G is computed (e.g., again on a classical computer) from the results 
of all measurements: V'(Co) • • • ,^k-i)- The probability measure A(f) is its 
distribution. For details, see Heinrich (2001). 

The error of A is defined as follows: Let < 9 < 1, f e F, and let C be 
any random variable with distribution A{f). Then put 

eiS,A,f,e) = inf {£ I P{||5(/) - C|| > s} < 9} . 

Associated with this we introduce 

e{S,A,F,e) = supe{S,A,f,e), 

e{S,A,f) = e{S,A,f,l/4), 
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and 

e{S,A,F)=e{S,A,FMA). 
The n-th minimal query error is defined for n G Nq as 

e'^{S,F) = mi{e{S,A,F) \ A is any quantum algorithm with nq{A) < n}. 

This is the minimal error which can be reached using at most n queries. The 
query complexity is defined for e > by 

compf(5,F) = 

mm.{nq{A) \ A is any quantum algorithm with e(5, A, F) < e}. 

The quantities en{S,F) and compf (S", are inverse to each other in the 
following sense: For all n G Nq and e > 0, en{S,F) < e if and only if 
compf J (5, -F) < n for all ei > e. Thus, determining the query complexity is 
equivalent to determining the n-th minimal error. Henceforth, we will deal 
only with en{S,F). 



2 The Main Result 

Let iV G N and set D = Z[0, N), K = K, G = K. For 1 < p < oo let L, 
denote the space of all functions / : Z) — > R, equipped with the norm 

/ N-l 

V i=0 

if p < oo and 

W/Wln = max \f{i)\. 

" 0<i<N--l ' ^ 

Define Sn ■ Lp — > R by 

N-l 

and let 

F = B^:={feL^\\\fhN<l}. 

Let us summarize the known results about the order of en{SN,Bp) (and 
thus the query complexity of computing the mean of p-summable sequences) 
in Theorem ||. The case p = oo is due to Grover (1998), Brassard, H0yer, 
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Mosca, and Tapp (2000) (upper bounds) and Nayak and Wu (1999) (lower 
bounds). The results in the case 1 < p < oo are due to Heinrich (2001). 
Note that throughout the paper we often use the same symbols for possibly 
different constants. Also, log always means log2. 

Theorem 1. Let 1 < p < oo. There are constants cq, 01,02,03 > such 
that for all n, iV G N with 2 <n< oiN , 

C2n-^ <eUSN,B^) <03n-^ if 2 < p < 00, 
C2n-i < e^(5iv,^^) < csn-^og^/^ nloglogn, 

and 

C2n-2(i-i/p) < eliSN.B^) < can^^Ci-i/p) -y. 1 < p < 2, n < coViV. 

The case 1 < p < 2, n > cq^/N was left open. We will settle it here by 
proving 

Theorem 2. Let 1 < p < 2. There are constants 00,01,02,03 > such that 
for all n,N £N with cqVN <n< oiN, 

C2n-2/PAr2/p-i < eiJSN,B^) < cgn-^/PiVS/p-i max(log(n/^), 1)2/p-i. 

It is interesting to mention the consequences for the case p = 1 sepa- 
rately: 

Corollary 1. There are constants 01,02,03 > such that 

C2<eUSN,B^) < 1 

ifO<n< a/]V, and 

02n-^N < el{SN,Bf) < csn'^N max{log{n/VN),l) 
ify/N <n< ciN . 

Hence the decay essentially starts only beyond \fN . Note that the cor- 
responding quantities for the classical deterministic and randomized setting 
remain ^2(1) also in the range ^/N < n < oiN, see Heinrich and Novak 
(2001). 

Combining this with the respective result in Theorem |l], we can cover 
the full range n < oiN. This result is a direct consequence of Theorems |^ 
and ^ and the monotonicity of en{SN, Bp) in n. 
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Corollary 2. Let 1 < p < 2. There are constants 01,02,03 > such that 
for all n, G N with n < ciN , 

C2min(n-2(i-Vp),„-2/PAr2/p-i) < eUSN,B^) 

< C3 min(n-2(i-i/p)^„-2/PAr2/p-i)max(log(n/ViV),l)2/P-i. 

The following two sections contain the proof of Theorem ^. 

3 Upper Bounds 

For any M G N we define 

SN,Mf = ^ Yl ■^(^) 

ieZ[0,Af),|/(i)l<M 

and ^ 

•S'tv.a// = SnJ - SN,Mf = ^ f{i)- 

ieZ[0,N),\f{i)\>M 

Proposition 1. Let 1 < p < 00. Then there is a constant c > such that 
for all n, M,N with 

n > cM-P/2^max(log(M-PiV), 1) 

we have 

Proof. It is easily verified that 

^n(^'n,M:Bp) = 

(we use the queries just classically to obtain the values of the f{i) up to 
any required precision and compute the sum classically). It follows that, 
modifying 0, if necessary, it suffices to prove the result for 

M > Mq. (6) 

We will specify Mq later on. Furthermore, we may also asssume that 

MP < N, (7) 
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because otherwise S'^ Mf — ^ f ^ ^oi^N m) — ^- 

m' = [log N] . (8) 

First we define a quantum algorithm from Bp to Z[0,2™ ) x R. To 
specify its quantum query, fix any m" > m' + 1 and define the mapping 
/3 : R ^ Z[0, 2"*") by setting for z e R 

2™"-^ if \z\ < M 

^2m"-m'-l(^ + 2"*')J if M < \z\ < 2"*' 

2™" - 1 if 2 > 2™' 

if 2 < -2"*'. 



It follows that for M < |z| < 2"^', 

_2m' ^ 2-"^"W+i^(^) < ^ < _2rn' + 2-"*"+"^'+i(/3(z) + 1), (9) 

and 

f3{z) = T""-^ if and only if \z\ < M. (10) 
In connection with this definition let us mention that for / € ^ 

< iV^/f < iV < 2"*' (i = 0,...,iV-l). (11) 

Put Z = Z[0,A^), let r : Z ^ Z[0,2"'') be the identical embedding, m = 
m! + m", and define the query by 

Q = {m, m', m" , Z, r, /3). 

Let = Hm' <8) -ffm" ) let 

(zeZ[0,2"*'), xGZ[0,2"*")) 

be the respective representation of basis states. Let Wq G U{H^i) be the 
Walsh-Hadamard transform, and let Xq G U{H^i) be defined by 



if i = Q 
\i) otherwise. 

Consider the following unitary transforms on H^, defined by: 

W\^\x) = {Wo\i))\x), 
X\i)\x) = iXo\i))\x), 

\i) \x) if ieZ andx^ T^"-^ 
\i) \x) otherwise, 
J |i) \x) = \i) \Qx) . 



T\i) \x) 
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Here Qx stands for (2™" — x)mod2'"". Note that Wq^ = Wq, and hence 
W-^ = W. For / G put 

= WXWQfJTQf. (12) 

Denote 

Df = {i\ieZ, \f{i)\>M]. 
It follows from the definitions above and from (p!oh that 



QfJTQfV)\^) = [ otherwise 

Aq will be an algorithm with one measurement. We define its unitary trans- 
form as 

QfYf^W, (13) 

where L S N will be specified later. The starting state will be |6o) = |0) |0), 
and the mapping : Z[0, 2™') x Z[0, 2"^") Z[0, 2™') x R will be given by 

(^(i, x) = {i, -2™' + 2-""+'"'+ix). (14) 

This completes the definition of algorithm ^o- Clearly, Yf is the Grover 
iterate for the set Df, and the whole algorithm is Grover 's search algorithm 
(Grover, 1996), or amplitude amplification, in the terminology of Brassard, 
H0yer, Mosca, and Tapp (2000), with respect to the Hm' component, fol- 
lowed by one more query Qf. Observe that by (^) and (11) each run of the 
algorithm Aq produces a pair {i,y) G Z[0,2"" ) X R with 

y</(i)<?/ + 2-'""+™'+i if teDf (15) 

and 

y = if and only if i < N and i ^ Df. (16) 

The final algorithm A is defined as ^(j4q ), which means that we repeat Aq 
L* times and compose the outputs by the mapping 

V^: (Z[0,2"') xR)^* ^R, 

see Heinrich (2001), Section 2, for a formal definition. The number L* € N 
will be specified later. The mapping ip is defined as follows: Let 



'-^ e (Z[0,2'"') X R) 
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be the outputs of the L* runs of Aq. We exclude all pairs with i£ ^ Df (which 
amounts to checking if i > or y = 0, by (|l6|)), as well as all repetitions 
of any i£ € Df (by a suitable sorting algorithm). For the remaining set we 
add the second components and divide by N (if the remaining set is empty, 
we output 0). 

Now we show that with a suitable choice of the parameters m", L, L* , the 
algorithm outputs S'j^ j^jf with error at most 2""* '^'^ +^ with probability at 
least 3/4. This follows from ( p!5| ) if we prove that with probability at least 
3/4 the set of remaining indices equals Df. If Df = 0, this is trivial, so we 
assume Df^%. First we analyze ^o- Denote = \Df\, hence /// > 1, and 
let < 61/ < 7r/2 be defined by 

sin^^/ = 2-"V/- (17) 

Finally, let 

i&Df 

and 

jeZ[0,2™')\D/ 

By the analysis of Brassard, H0yer, Mosca, and Tapp (2000), relation (8), 

y/^y |0) |0) = (2--V/)-'/' sin((2L + 1)6 f) |^/,i) |0) + 

(1 - 2— cos((2L + 1)9 f) |^/,o) |0) , 

(where the second term is replaced by if ^/ = 2™' ) . It follows that for any 
io S Df, the algorithm Aq outputs (io, Pifi^o))) with probability 

gi,=fifsm\{2L + l)9f). (18) 

In the sequel we use the elementary relation 

2x/7r < sinx < X (x G [0, 7r/2]). (19) 

Since / G Bp , we have 

N-'^MP\Df \ < 1, 

hence 

^if = \Df\< M'PN (20) 
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and 

2-""'^/ < M-PiV2-™' < M-P. 
Therefore, by (0) and (|l^ 



and hence 



Now we put 



and define L by 



< (211 



Mo = [e^/P] (22) 



L=[3-^MP/2j. (23) 
Since we assumed M > Mq, we get from (p^) and (p3|) , 



1 < Imp/^ < L < (24) 
6 3 

It follows from (|2l|) and (||) that 

(2L + 1)6*/ < 3L(9/ < 7r/2. (25) 
On the other hand, by (El) and (O), 



(2L + 1)^/ > 2L% > -MP/2sin0j = -Mp/2(2-™'^^)1/2_ 
3 3 



From (^), (|l8|), ( ^5|) and the relation above, 

- 97r2 

> -^MPN^^ = C2MPN-^, 

- 97r2 

where in the last line we used and set C2 = 2/(9'ir'^). It follows that after 
L* repetitions of algorithm Aq the probability of {io, (3{f{io))) not being 
among the results is 

< (1 - C2MfiV-i)^* < e--2MPiV-iL*^ 
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where we used that 1 + x < for x G R. The probabihty that at least one 
io £ Df is not among the results is 

where we used (^0|). Now we choose L* in such a way that this probability 
is not greater than 1/4. This requires (recall that log means log2) 

(C2 loge) M'PN-^L* > \og{M-PN) + 2, 

which is satisfied if 

3 



C2 log e 



M-PiVmax(log(M'PAf), 1) 



We put C3 = 3/(c2loge) and observe that the above combined with (^) 
implies 

L* < (C3 + l)M-PiVmax(log(M-PiV), 1). 
Together with (^), this implies that algorithm A makes 

(2L + l)L* < 3LL* < (c3 + l)M-P^'^Nmax{log{M-PN), 1) 

queries to compute S'^ j^^f up to error 2"™ with probability at least 

3/4. Since m" was arbitrary, the result follows. □ 

We need to express M in terms of n and A'^: 



Corollary 3. Let 1 < p < oo. There is a constant c > 1 such that for all 
n, M, N eN, 

whenever 

M > c{N/nf/P max(log(n/\/]V), if/^. 
Proof. Let cq be the constant from Proposition |l|. We put 

c = max((2co)^/P,l). (26) 

Assume 

M > c(iV/n)2/Pmax(log(n/\/iV),l)2/P. 

It follows that 

M-P/'^N < c^P/2^/max(log(n/\/iV), 1). (27) 
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Squaring and dividing by N gives 

M-PN < c"ViV"Vmax(log(n/\/]V), 1)2, 

and hence 

max(log(M-PiV),l) 

< max log(c~P) + 2 \og{n/VN) - 2 log ( max(log(n/\/]V), 1)) , 1^ 

< 2max(log(n/\/iV),l). (28) 

(H, (|8D and dH) give 

coM-P/^Nmax(log{M~PN),l) < IcQc'^l'^n < n, 
which, by Proposition ||, impUes 

□ 

Proposition 2. Let 1 < p < 2. There is a constant c > such that for all 
k,n,N e N, 

4(S^,2^,0 < c(2(i-P/2)'=n-i +2'=n-2). 

Proof. This is a direct consequence of the method of proof of Theorem 1 
in Heinrich (2001). For the sake of completeness, we recall some key steps. 
Since trivially en{Spf 2''T^p) ^ 1 for ^-H ri G Nq (just use the zero algorithm), 
it suffices to prove the result under the assumption 

n > 2(i-P/2)^ (29) 
Define 5^*^ : — > R for £ = 0, . . . , fc, o" = 0, 1 as 

S'jff = (-l)-2-^iV-i 

2^-l<(-l)'^/{i)<2< 

if £ > 1 and 

o<(-i)'^/W<i 

It is shown in Heinrich (2001) (based on the counting algorithm of Brassard, 
H0yer, Mosca, and Tapp, 2000), that there is a constant c > such that 



13 



for each choice of I'e, n£ £ N (i = 0, . . . ,k), there are algorithms Ai^a- = 
0, . . . , fc, o" = 0, 1) with ng(yl^^o-) < i^en-e and 

(use the relation following (27) in Heinrich, 2001, together with (21) and 
(22) of that paper). Now choose 



and 
Due to (12 



2-{l/2-p/4){fc-£)^ 

U£ = [21og(A;-£ + l)] +4. 



Let the algorithm A be defined by 

0<£<fe 

CT = 0,1 

(We refer again to Heinrich, 2001, Section 2, for a formal definition.) Taking 
into account (|3^), it follows that 

k 



n,{A)<2Y,mogik -£ + !)]+ A) 

1=0 

Moreover, since 



-(l/2-p/4)(fc-^) 



n 



< cin. (31) 



£=0 



£=0 



we get 



- ^E (2(^~P/2)^+(^/2~P/^^(^"^)n"^ + 2- 



/2-p/4){fc-f)„-l ^ ^f+(l-p/2)(fc-^)^-2 



^=0 



f=0 



- ^E (2(^/2-p/4)(fc+<?)^-i _^ 2^-p(*=-^)/2n-2 



< C2 2 



which together with (^) and a suitable scaling of n implies the desired 
result. □ 
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Theorem 3. Let 1 <p <2. There are constants cq,c> such that for all 
n,N £N with n > cqVN 

el{SN,B^) < cn-2/Piv2/p-imax(log(n/^/iV),l)2/f-i. 
Proof. First note that 

e%{SN,B^) = 0. (32) 

Next observe that it follows readily from Lemma 3 in Heinrich (2001) (re- 
ducing the error probability by repeating the algorithm and computing the 
median) that there is a constant cq G N such that for all n,k,N £ N, 

e?„J5;v,0 < e^^{S^^,,,B^)+eUS'^,„B^). (33) 

Now let n satisfy 

< n < (34) 
and choose A: G N in such a way that 

2'=^^ < ci(iV/n)2/Pmax(log(n/^), if/^ < 2^ 
where ci > 1 is the constant from Corollary |3|. Consequently, we have 

Moreover, with C2 being the constant from Proposition |2|, 

(^{Sjq2k,Bp) 

< C2(2(i-P/2)'=n~i + 2'=n-2) 

n \ p^^ p/2) 



< C3 \^{N/n)¥'^-''^^^n~'^ max (^log 1^ 

log--^,lJ 

= C3 ^iV2/P-in-2/P max (^log l) ' 

+iV2/p^-2/p-2 (log 1^ . (36) 



2/pN 



2/p-l 
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Using (again) x > ln(l + x) for x > —1, we have 

r? ^ ( r? n 2, n , n 
> In hi > 2 In ^= = log ^= > log 



N \N J ViV loge ^/N VN 

Consequently, recalling our assumption n > \/]V, we get 



2 

n / , n 



->max(log^,l 



and therefore 

/ \ 2/p-l / ^ 2/p 

j^2/p-i^~2/p max ( log 1 > Ar2/p„-2/p-2 ( log 1 

From (H), (il), (|36|) , and the relation above we get 

< 2c37V2/p-l^-2/p J^^Qg 

for all n with \/iV < n < N. With a suitable scaling of n, the result follows 



2/p-l 

(37) 



from (HD and (H). □ 



4 Lower Bounds 

We need some general results from Section 4 of Heinrich (2001). Let D and 
K be nonempty sets, let L G N, and let to each u = {uq, ■ ■ ■ , ul-i) G {0, 1}^ 
an /u G J^{D, K) be assigned such that the following is satisfied: 

Condition (I): For each t £ D there is an ^, < ^ < L — 1, such that 
fu{t) depends only on U£, in other words, for u, u' E {0, 1}^, = ti^ implies 

Ut) = fu'{t). 

Define the function Q{L,e,£') for L E N, < £ / f < L by 



The following was proved in Heinrich (2001), using the polynomial method 
of Beals, Buhrman, Cleve, and Mosca (1998) and based on a result of Nayak 
and Wu (1999): 
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Lemma 1. There is a constant cq > such that the following holds: Let 
D, K he nonempty sets, let F C J^{D, K) he a set of functions, G a normed 
space, S : F ^ G a function, and L € N. Suppose {fu)u&{o,i}'- ^ ^{D,K) 
is a system of functions satisfying condition (I). Let finally < i ^ £' < L 
and assume that 

fu^F whenever \u\e {£,£'}. (39) 

Then 

eUS,F) > imin{||5(/„) | H = £, \u'\ = £'} (40) 

for all n with 

n<cog{LJ,i'). (41) 

The next result contains lower bounds matching the upper ones from 
Theorem ^ up to a logarithmic factor. 

Theorem 4. Let 1 < p < 2. Then there are constants co,ci,C2 > such 
that for all n, G N with cqVN <n< ciN , 

Proof. Let cq be the constant from Lemma |l|, and let 

ci = co/Vu. (42) 

By assumption, 

coVn <n< ciN. (43) 

We set 

L = N, e= [2co ^n^AT-i] , / = ^ + 1. (44) 
It follows from (El) that £ > 2. Moreover, from 



n < CQ^frnfi (45) 
and, taking into account that ^ > 2, 

ll2<l-\< 2co2n2iV-\ 
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hence, by (|4|) and (|3|), 

£ + 1 < 3^/2 < Qc^'^n^N-^ < GcQ^cfN = N/2. (46) 
We have, by (H), (||) and (|4|), 



n<coVeN/2<co min Vj(iV - j) < co^?(L, ^, ^')- (47) 



Now we define ipj G Lp (j = 0, . . . , L — 1) as 



We have 



{£ + l)-^/PN^/P if i = j 
otherwise. 



SniPj = {e + i)-^/'PN^/p^\ 

For each u = {uq, . . . , ul-i) G {0, 1}^ define 

L-l 

fu = ^Uj^j. (48) 
i=o 

Since the functions ilij have disjoint supports, the system {fu)u£{o,i}^ sat- 
isfies condition (I). Moreover, fu S whenever \u\ = £,£ + 1. Lemma ||, 
relation (47) and the left and middle part of ( p^ ) give 

eUSN,13^) > ^mm{\SNfu-SNfu'\\\u\=£,\u'\=£ + l} 

= i(^ + lyi/p^i/p-i > 1 (g^-2^2^-i)-i/p^i/p-i 



«— n~2/PAr2/P-i. 
2 • 6Vp 



□ 



5 Comments 

Let us first mention that there remains another gap in the order of the 
quantity en{S]\fiBp) in all the results of Theorems |l|, ^, and Corollaries ||, 
namely, the region ciN < n < N . As we mentioned before, we have 
en{SN,Bp) = for n> N (classical computation of the sum). Hence filling 
this gap means determining how fast ekiS^ ,3^) goes to zero in the region 
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close to classical computation. We did not consider this problem further. 
It is theoretically interesting, but one should also mention that its solu- 
tion would not say much about the speed-up due to quantum computation: 
With an effort, just by a constant factor higher, the problem can be solved 
with the same error (in fact, even up to any needed precision) by classical 
computation. 

Finally, we discuss the cost of our algorithm in the bit model of computa- 
tion. Here we assume that both N and n are powers of two. The algorithm 
behind Proposition || and Corollary ^ needs 0{nm") quantum gates (see 
Nielsen and Chuang, 2000, Chapter 4, for basics on quantum gates), 0{m") 
qubits, and makes 0(n^A^~^/ max(log(n/\/]V), 1)) measurements to reach 
error 0(2'°g^~™' ). The bit cost of the classical computations is negligi- 
ble as compared to the number of quantum gates: We need 0{n?N^^m") 
classical bit operations to sort out the wrong elements and to add the right 
ones. The bit cost of the algorithm in connection with Proposition |2| was al- 
ready analyzed in Heinrich (2001). It amounts to 0{nlogN) quantum gates, 
0{logN) qubits, and 0{k\ogk) (which is ©(log n log log re)) measurements. 
The number of classical bit operations is O (log re log log n log A^), and thus, 
again dominated by the number of quantum gates. Summarizing this for the 
algorithm of Theorem we see that we can implement it with 0(relog A^) 
quantum gates, on 0{logN) qubits, and with 

©(n^Af" V max(log(re/\/iV), 1) + log(iV/re) log log(7V/re)) 

measurements. Thus the quantum bit cost differs by at most a logarithmic 
factor from the quantum query complexity. 
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